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1. INTRODUCTION 

A signal can be completely rebuilt using a linear reconstruction method if it is sampled at equal intervals 
and bandlimited by less than the sampling rate. The frequency range and sample rate spanned by the signal are 
used to reconstruct the signal. The signals can be captured and modeled by several mathematical approaches. One 
such kind of approach is the compressed sensing (CS) technique [1]-[8]. CS assumes the signal to be sparse. For 
a two-dimensional signal, most of the entries in its vector space are zeros. So, it is called as sparse vector space 
for limited dimensional signal. Sparse vector doesn’t have itself to be consist of signal of interest but it can be 
equally well represented by some basic functions like wavelet or fourier transforms. 

Let’s consider a set of vectors {x} that are linearly dependent and if there are scalars that not at all 
zeros is given by (1) X; Aixi = 0. Then set of vectors {xi} are linearly independent and orthogonal. If the inner 
product of any two vectors in the set is zero, then it is called orthogonal. Sparsity is an important and more 
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powerful concept in the field of mathematics and computer science. Consider an N-dimensional x, it can be 
expressed by a linear super position of K<<N elementary signal quantities for some error e and e(k) called 
sparse for the vector a. K quantize yw, are collected from a dictionary in turn which is represented by N x D 
matrix with D>=N [2]. 


x = Dear CHW in ||x — Xgll2 < e(k) (1) 


The exact sparsity can be represented by y=¥x, [9]-[11] where ¥ is an orthonormal basis, x is the input 
signal. Some of the dictionary representations are fast fourier transforms (FFT), discrete cosine transform (DCT), 
and discrete wavelet transform (DWT). FFT [12] of the signal xx of the signal yx is given by (2). 


X= U/sqrt(n) Die Yk exp(—2n(j — 1)(k — 1)i/N) (2) 


By assuming the sparse version of signal X [13] in which all significant entries are set to zero. We can use this 
sparse vector Z to approximate the original signal. 6Y =Y—Y = F(X —Z)we can make of Z to 
approximately predict y. Z > Y ~ Y > y. As we know the locations of non-vanishing entries in the sparse 
vector, only fractions of samples in Y are used to recreate Z [14] as in (3). 


Then y™ = FmxmZ™ = FinxnZ (3) 


Where y™ € R™,Z™ € R™, F nxm ermXm Fy yeRMXN, are got by removing the rows/columns of the vanishing 
entries in Z. And so, we can recreate by sampling m<N points at specific locations which may go faraway less 


than the Nyquist rate [15]. y" = FmxnZ, in which the sample in Z are less than the length of the sparse vector. 
The relationship between the discrete representation of y and fourier transform is given by X=FY where F 
being the DFT dictionary, X = {X,}W_, E€ R” and Y = {Yp}; E€ R”. The inverse way of producing the signal 
Y = FX, where F is transpose of the F (inverse of DFT) and X is called the compressible vector or sparse 
vector. In this case of signal X, the signal energy is not distributed equally over the spectrum. So, it is 
approximated by reduced representation. DCT [16] is considered as the discrete-time version of the fourier 
cosine series as shown mathematically (4). DCT estimates the real approximation of the signal with less 
coefficients. 


Dj- CNINCH),CH ES TN" pe, y) cos VY, cos PtP] (4) 


Otherwise, C(u)=C(v)=1. The pixel in the image represented by the M x N matrix is p (x, y). DCT is 
applied on a block of size N, and f(x, y) is the N x N input image. The encoded or transformed image of f(x,y) is 
F(u,v). DWT [17], [18] decomposes the signal into approximation and detailed subbands. The detailed 
subbands are having the much information about the image. They are called as high significant components of 
the signal. Approximate subbands are the least significant components of the mages. The least significant 
components are nothing but the edges present in the image. There are different types of wavelet families like. 
Haar, Dabhache, and symmetric. We have experimented with Haar wavelet transforms. The general wavelet 
transform is (5) and (6). 


Fuign pass(™, n= fS f Venn (Eat (5) 


Frow pass (m N)=J_.. f (xX) Pin ny tdt (6) 


High pass filter Fyign pass(™,n) passes through high-frequency components and low pass filter 
Flow pass(M,N) passes through low-frequency components. * is a conjugate symbol. In (5) represents the 
wavelet function given by Y’ and in (6) represents the scaling function given by ®’. f(x) is the input signal. 
The input signal is represented by f(x). The wavelet transform has some advantages: i) wavelet transforms at 
higher levels avoid blocking artifacts, ii) well matched to Herpes simplex virus (HSV) characteristics, and iii) 
wavelet provides higher level of decomposition. 

Let A denote the m x N sampling matrix with the restricted isometry criterion (RIC) 62s <=c [19]. Let 
consider a signal f that is sparse, x has only K non-zero entries. f=y x, where w is some transform domain [20], 
[21]. The measurement vector can be obtained as y = ġf = dx = Ax . A is called the measurement matrix 
and @ is called the sampling matrix. The sparse vector can be reconstructed by solving the linear equation 
y=Ax [22]-[24]. 
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2. PROPOSED METHODS 

The compressive sampling matching pursuit (CoSaMP) [25], [26] is a compressed sensing technique 
based on orthogonal matching pursuit. It must fulfill the following requirements: i) it must begin with the least 
sample data possible. ii) It should progress by taking into account samples from various sampling strategies. 
iii) It must accept all samples that have been combined with noise and must be resilient. Novel CS techniques 
DCTCoSaMP, FFTCoSaMP, and DWTCoSaMP are proposed in this paper to recover the original signal from 
the sparse or incomplete signal components. Finally, using the peak signal-to-noise ratio (PSNR), structured 
similarity indexing method (SSIM) and root mean square error (RMSE) approaches, the restoration of magnetic 
resonance imaging (MRI) images is analyzed and compared. 


2.1. Orthonormal bases 

Assume set B. B has a set of vectors, viz. B=v1, v2, ......vk, each of which has a length of 1. ||vil|=1 
for i=1, 2,...k. or vi. vi=1 or ||vil]2=1for i=1,2,3......k. All of the vectors in B have a length of one, indicating 
that they all have been scaled and are unit vectors. Any two vectors in the collection have a dot product of zero, 
i.e. vi. vj=0 for i and j. These vectors are all orthogonal and have the same length and are normalized. As a 
result, set B is considered too as orthonormal. The term orthogonal points that all of the set's components are 
orthogonal and normalized. 

CoSaMP is an orthogonal matching pursuit-based (OMP-based) encoding/decoding approach. 
Because of one of the compressed sensing techniques, CoSaMP has gotten a lot of interest today. The CoSaMP 
method accepts a K-sparse signal array as input. The measurement grid should be accepted. Finally, it should 
also accept the criterion for halting. The CoSaMP [27]-[29] is a greedy matching pursuit algorithm. 


2.2. Algorithm CoSaMP (M, phi, s) 

The algorithm begins with a basic signal approximation, in which the starting leftover signal equals 
the unknown target signal. It functions through 5 important phases throughout each iterative process: 
i) identification: from the available samples, the algorithm creates a substitute of the residual and locates the 
significant components of the substitute; ii) support merger: the newly discovered components are combined 
with the components found in the present estimation; iii) estimation: to estimate the target signal on the merged 
set of components, the method performs a least-squares problem; iv) pruning: by preserving just the biggest 
elements in this least-squares signal approximation, the method creates a new approximation; 
v) sample update: finally, the samples are updated to represent the residual signal, which is the portion of the 
signal that has not been estimated or approximated. 

Figure 1 depicts the recommended architecture for the proposed approaches. The system receives an 
MRI image as input, and uses dictionary extraction methods including FFT, DCT, and DWT to leverage 
sparsity factors. The low and high-frequency components are included in the dictionary. Only low-frequency 
components, which signify incomplete or sparse components, are retained. Novel CS techniques DCTCoSaMP, 
FFTCoSaMP and DWTCoSaMP are proposed in this paper to recover the original signal from the sparse or 
incomplete signal components. Finally, using the PSNR, SSIM, and RMSE approaches, the restoration of MRI 
images is analyzed and compared. 


Imput: Significant coefficients to be selected M, random measurement matrix Phi, 

sparsity level s. 

Ouput: The target signal's s-sparse approximation 

Initialisation: fo €0 // initial approximation, yoy //Current samples = input 
samples, t€0 


loop 
1. te tril // Increase the number of cycles each time. 
2 5 € Mty® // Create a list of intermediates. 
3. B €© supp (62s) // Make a comprehensive collection of substitutes. 
4. BEB U supp(f*?) // Promote a merger 
5. PIBE € (Mrt)* Bt. // Least-squares signal estimate 
6. PIB: €0 
7 ft € b,x // Shear/trim in order to get the next estimate 


8. y? € y - Aft // Renew the samples: 
till halting criterion holds 
fe £ 


The block diagram for DCTCoSaMP/FFTCoSaMP/DWTCoSaMP as shwon in Figure 2. The CoSaMP 
algorithm takes the input iage, where it calculates the automatic correlation, it means, it selects only the sparsity 
features present in the signal. And it selects the current set of preselected atoms (only a few k non-zero 
coefficients). Then merge sparse vector with measurement matrix to form the candidate atom set. Then least 
square methods are used to get the signal residuals. Shear or prune to obtain the supporting atom set in order to 
get the better estimation. Update the current final estimate and calculate the current residuals until stopping 
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criteria is met. If stopping criteria is not met, update the current residual and support set. And the whole process 
continues further. 


Under sampled 


Input MRI Image K Sparse Data 


Sparse Dictionaries obtained 
via Transform Domains 


Reconstructed Image is 
Verified via PSNR, SSIM Apply zero to the 
and RMSE Place of Sparse Data 
Position 


Novel Techniques 
DCTCoSaMP, 
- 


Figure |. Architecture for the proposed methods 


Figure 2. Block diagram for FFTCoSaMP/DCTCoSaMP/DWTCoSaMP 


3. EXPERIMENTAL RESULTS 

In this article, six different types of MRI images like Glioma Tumor images, Brain Tumor images, 
Brain Tumorless images, no Brain Tumor images, Non Cancerous images, and Cancerous images as shown in 
Figure 3. The experimentation is conducted on about 2,000 MRI images. The dataset has been taken from some 
websites like www.kaggle.com and the state-run Victoria hospital, Bengaluru. We have experimented on 
different dimensions of measurement matrices like 256x210, 256x180, and 256x155. Among all these 
dimensions of measurement matrices, the dimension 256x210 is exploiting the best reconstruction and 
acceptable measurement matrix. For the 256x210 measurement matrix, the comparative results for different CS 
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techniques like FFTCoSaMP, DCTCoSaMP, and DWTCoSaMP are tabulated and analyzed through 
experimentation methods. 

Table 1 compares PSNR, SSIM, and RMSE values obtained from thorax MRI images and compared 
to some research articles using various compressed sensing approaches. The proposed methods FFTCoSaMP, 
DCTCoSaMP, and DWTCoSaMP provide better PSNR, SSIM, and MSE values than all of the references 
found in the literature review. The results of experimentation are tabulated. 

CS technique has been discussed in various papers as per the references mentioned in [26], [30]-[33]. 
Where PSNR values for different CS techniques including CoSaMP algorithm are discussed as per the 
experiments from various researchers, it is observed that PSNR values for my proposed methods are higher 
than the other CS techniques of the literature survey. Similarly, SSIM values for different CS techniques 
including CoSaMP algorithm are discussed as per the experiments from various researchers. It is observed that 
SSIM values for my proposed methods are higher than the other CS techniques of the references [34]-[37]. 
Finally, RMSE values for different CS techniques, are discussed as per the experiments from various 
researchers. It is observed that RMSE values for my proposed methods are lower than the other CS techniques 
of the literature survey [34]-[38]. 


Table 1. PSNR, SSIM, and RMSE values for the 


proposed methods 
ieee eae Oe 
lioma T j j s 
G eae Brain Tumor MRI Brain a PSNR FFTCoSaMP TA 
= tine Ova. » 
A à - DCTCoSaMP 36.33 
í a Fe DWTCoSaMP 37.16 (DB4) 
i j 38.12 (Coif3) 
5 38.5 (Sym8) 
DaS a 
Meningioma Cancerous MRI Non-Cancerous SSIM FFTCoSaMP 0.7 
Tumor MRI MRI DCTCoSaMP 0.73 
DWTCoSaMP 0.81 
y : be ot RMSE FFTCoSaMP 0.66 
Figure 3. MRI images used for exploitation of DCTCoSaMP 0:53 
various sparse dictionaries DWTCoSaMP 041 


Table 2 shows the comparison of different compressive sampling techniques w.r.t. PSNR, SSIM and 
RMSE for different MR images like Glioma Tumor images, Brain Tumor images, Brain Tumorless images, no 
Brain Tumor images, non Cancerous images, and Cancerous images. As mentioned, the experiment is 
conducted on 2,000 MRI of of different categories of Glioma-Tumor images, Brain Tumor images, Brain- 
Tumorless images, no_Brain_Tumor_images, non-Cancerous images, and Cancerous MR images. The average 
value of PSNR for all these images w.r.t. CS technique DWTCoSaMP is 33.0347, which is high when 
compared to other CS techniques like FFTDCoSaMP and DWTCoSaMP with PSNR values 26.6229 and 
27.7749 respectively. Similarly, the average value of SSIM for all these images w.r.t. CS technique 
DWTCoSaMP is 0.6373, which is almost equivalent to SSIM value 0.7136 of FFTCoSaMP and 0.4506 of 
DCTCoSaMP. Lastly the average value of RMSE for all these images w.r.t. CS technique DWTCoSaMP is 
0.0019, which is less when compared to other CS techniques like FFTCoSaMP and DCTCoSaMP with PSNR 
values 0.0451 and 0.0036 respectively. 


Table 2. Comparison of different compressive sampling techniques w.r.t. PSNR, SSIM, and RMSE for 
measurement matrix of dimension 210 X 256 


s Metrics PSNR SSIM RMSE 
N MRI images FFTCo DCTCo DWTCo FFTCo DCTCo DWTCo FFTCo DCTCo DWTCo 
g SaMP SaMP SaMP SaMP SaMP SaMP SaMP SaMP SaMP 

1 Glioma Tumor 26.7953 27.9552 32.9565 0.7021 0.4644 0.5938 0.0016 0.0019 0.0006 
images 

2 Brain Tumor 26.1371 27.3214 32.7034 0.8534 0.4914 0.8423 0.2190 0.0021 0.0006 
images 

3 Brain Tumorless 26.7877 27.7369 33.2316 0.4284 0.4251 0.5759 0.0016 0.0019 0.0006 
images 

4 No Brain Tumor 26.9530 28.2062 33.3209 0.7264 0.4311 0.5891 0.0015 0.0053 0.0069 
images 

5 Non-Cancerous 26.4417 27.6548 32.9611 0.8575 0.4412 0.5855 0.0017 0.0068 0.0006 
images 

6 Cancerous 26.0886 27.3342 32.6845 0.4940 0.4915 0.8040 0.0018 0.0021 0.0049 
images 
Average 26.6229 27.7749 33.0347 0.7136 0.4506 0.6373 0.0451 0.0036 0.0019 
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As tabulated in the Table 3, DWTCoSaMP is the best when compared to other CS techniques like 
FFTCoSaMP and DCTCoSaMP. The average values of PSNR for DWTCoSaMP are 32.9763, 32.1498, and 
28.7848 for different measurement matrices like 210 x 256, 180 x 256, and 155 x 256 respectively. Similarly 
The average values of PSNR for DCTCoSaMP are 27.7014 dB, 26.1703 dB, and 24.2075 dB for different 
measurement matrices like 210 x 256, 180 x 256, and 155 x 256 respectively. 


Table 3. PSNR based comparison of different compressive sampling algorithms under different dimensions 
of measurement matrices 


C> DWTCoSaMP DCTCoSaMP FFTCoSaMP 
algorithms 
SN eae 155x 180x 210x 155x 180x 210x 155x 180x 210x 
ai AX 256 256 256 256 256 256 256 256 256 
MRI images 
1 Glioma Tumor 29.4521 32.7025 32.9565 24.0214 26.4812 27.9552 29.4521 25.0786 26.7953 
images 
2 Brain Tumor 29.0123 31.8925 32.7034 24.1245 26.7891 27.3214 28.0123 24.9214 26.1371 
images 
P 3 Brain 29.4358 32.1245 33.2316 23.9754 26.2471 27.7369 29.4358 24.1234 26.7877 
s Tumorless 
N images 
R 4 No Brain 28.1245 32.1246 33.3209 23.8974 25.8974 28.2062 28.0245 25.3124 26.9530 
Tumor images 
5 Non-cancerous 27.8975 32.3791 32.9611 24.5687 25.4589 27.6548 27.8975 24.5698 26.4417 
images 
6 Cancerous 28.7864 31.6754 32.6845 24.6573 26.1478 27.3342 28.7864 25.1234 26.0886 
images 
Average PSNR 28.7848 32.1498 32.9763 24.2075 26.1703 27.7014 28.6014 24.8548 26.5339 


Table 3 shows the PSNR based comparison of different compressive sampling algorithms under 
different sizes of measurement matrices like 155x256, 180x256, and 210x256. As well, the average values of 
PSNR for FFTCoSaMP are 27.7014 dB, 26.1703 dB, and 24.2075 dB for different measurement matrices like 
210x256, 180x256, and 155x256 respectively. Similarly average values of SSIM for the same differnet 
measurement marices is tabulated and compared in Table 4. Likewise, average values of RMSE for the same 
different set of measurement matrices are tabulated and compared in Table 5. 


Table 4. SSIM based comparison of different compressive sampling algorithms under different dimensions of 
measurement matrices 


cS DWTCoSaMP DCTCoSaMP FFTCoSaMP 
Algorithms 
EN ee 155x 180x 210x 155x 180x 210x 155x 180x 210x 
nd ae 256 256 256 256 256 256 256 256 256 
MRI images 
1 Glioma Tumor 0.5645 0.5332 0.5938 0.4674 0.4744 0.4644 0.4650 0.5144 0.7021 
Images 
2 Brain Tumor 0.4919 0.7424 0.8423 0.4824 03914 04914 04754 05914 0.8534 
Images 
S 3 Brain Tumorless 0.4268 0.5453 0.5759 0.5251 0.5251 0.4251 0.5721 0.4251 0.4284 
S Images 
I 4 No Brain Tumor 0.4611 0.5861 0.5891 0.4411 0.3911 0.4311 0.5411 0.5311 0.7264 
M Images 
5 Non-Cancerous 0.4412 04855 0.5855 0.4342 0.4121 04412 0.4412 0.4412 0.8575 
Images 
6 Cancerous 0.4915 0.7940 0.8040 0.5115 04115 04915 0.4250 0.4915 0.4940 
Images 
Average SSIM 0.4795 0.6144 0.6651 0.4769 0.4343 0.4574 0.4866 0.4991 0.6770 


4. DISCUSSION 

About 2000 MRI images like Glioma Tumor images, Dian Tumor images, Tumor-Less images, Cancerous 
images and Non Cancerous images are used for experimentation. For all these images, average PSNR values for 
DWTCoSaMP, DCTCoSaMP, FFTCoSaMP are 32.9763, and 27.7014, 26.5339 respectively for the measurement 
matrix of 210x256. Similarly, SSIM values for DWTCoSaMP, DCTCoSaMP, FFTCoSaMP are 0.6651, 0.4574, and 
0.6770 respectively for measurement matrix 210x256. Lastly, RMSE values for DWTCoSaMP, DCTCoSaMP, and 
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FFTCoSaMP is 0.0017, 0.0034, and 0.0379 respectively for measurement matrix 210x256. In all data analysis made 
in different measurement matrices of 155x256, 180x256, and 210x255, CS technique DWTCoSaMP is the best 
compered to DCTCoSaMP and FFTCoSaMP and other CS techniques. 

Figure 4 shows the bar chart for the various PSNR values of Table 3. Figure 5 shows the bar chart for the 
various SSIM values of Table 4. Figure 6 shows the bar chart for the various RMSE values of Table 5. The charts 
show that for all different measurement matrices, the DWTCoSaMP provides a best result when compared to other 
CS techniques of the letrature survey as refered in various references. If the SSIM value is appearing towards 1 (one) 
or higher the PSNR values towards 100 dB or lower the RMSE values towards 0 (zero), there is a very good and 
better improvement in the reconstruction of MRI with respect to sparsity and incomplete signal. 


35,0000 
%—=r]- II - PSNR 
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SESS] SEAE S=E%||- J) SEAE EE SEAE SEA) SEAE 
20,0000 FEE EA SSF SEESE] AREA- 
e—a A iai lie zsif- a Zezi- Nszvi-| Azs- Ajj- =| (Bi 
15,0000 -Q=54}-|Q=24] =| EEA- -=+ 1- N= =| SEA Ale = 7 
yzl- ket] - kasije J ==#7 =A E jedmi lie al adi | pe 
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œ Glioma-Tumor Images = Brain Tumor Images = Brain-Tumorless Images 
#No_Brain_Tumor_Images I! Non-Cancerous Images - Cancerous MR Images 


II Average_PSNR 


Figure 4. PSNR based bar chart for the analysis of FFTCoSaMP, DCTCoSaMP, and DWTCoSaMP 
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155 x 256180 x 256210 X 256155 x 256 180 x 256210 X 256155 x 256 180 x 256210 x 256 
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Figure 5. SSIM based bar chart for the analysis of FFTCoSaMP, DCTCoSaMP, and DWTCoSaMP 
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Figure 6. RMSE based bar chart for the analysis of FFTCoSaMP, DCTCoSaMP and DWTCoSaMP 


Table 5. RMSE based comparison of different compressive sampling algorithms under different dimensions 
of measurement matrices 


cs DWTCoSaMP DCTCoSaMP FFTCoSaMP 
Algorithms 
SIN A 155x 180x 210x 155x 180x 210x 155x 180x 210x 
i ATX OSG 256 256 256 256 256 256 256 256 
MRI images 
1 Glioma Tumor 0.0329 0.0403 0.0006 0.0021 0.0036 0.0019 0.2004 0.1238 0.0016 
images 
2 Brain Tumor 0.0330 0.0223 0.0006 0.0224 0.0009 0.0021 0.1928 0.1237 0.2190 
images 
3 Brain 0.0332 0.0021 0.0006 0.0192 0.0029 0.0019 0.1271 0.0049 0.0016 
R Tumorless 
M images 
S 4 No Brain 0.0330 0.0013 0.0069 0.0539 0.0011 0.0053 0.1245 0.0004 0.0015 
E Tumor images 
5 Non- 0.0430 0.0005 0.0006 0.0210 0.0407 0.0068 0.0298 0.0005 0.0017 
Cancerous 
images 
6 Cancerous 0.0236 0.0011 0.0006 0.0234 0.0005 0.0021 0.1314 0.0011 0.0018 
images 
Average RMSE 0.0331 0.0113 0.0017 0.0237 0.0083 0.0034 0.134 0.0424 0.0379 


5. CONCLUSION 

In this research work for exploitation of sparsity under different compressed sensing techniques, 
DWTCoSaMP reveals the best method when compared to DCTCoSaMP and FFTCoSaMP and other traditional 
methods. By varying the dimensionalities of measurement matrices, DWTCoSaMP gives best results when 
compared to DCTCoSaMP, FFTCoSaMP, and other traditional CS techniques. For my future study, we would 
like to work on audio and speech signal type of data. So, we would like to work on these datasets. 
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